General formula for the four-quark condensate and vacuum 

factorization assumption 

Hong-shi Zong 1 ' 2,3 , Deng-ke He 1 , Feng-yao Hou 1 , and Wei-min Sun 1,2 
1 Department of Physics, Nanjing University, Nanjing 210093, China 
2 Joint Center for Particle, Nuclear Physics and Cosmology, Nanjing 210093, China and 
3 CCAST(World Laboratory), P.O. Box 8730, Beijing 100080, China 

Abstract 

By differentiating the dressed quark propagator with respect to a variable background field, 
the linear response of the dressed quark propagator in the presence of the background field can 
be obtained. From this general method, using the vector background field as an illustration, 
we derive a general formula for the four-quark condensate (0| : 9(0)7 At g(0)g(0)7 M g i (0) : |0). This 
formula contains the corresponding fully dressed vector vertex and it is shown that factorization 
for (0| : 9(0)7^(7(0)5(0)7^(7(0) : |0) holds only when the dressed vertex is taken to be the bare one. 
This property also holds for all other type of four-quark condensate. 
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It is well-known that one of the central problems of strong-interaction physics is in under- 
standing the structure of the ground state of quantum Quantum Chromodynamics, the QCD 
vacuum. One way to characterize this structure is by means of various vacuum condensates, 
such as the two-quark condensate, the gluon condensate and the four-quark condensate, etc. 
These condensates are essential for describing the strong interaction physics using the QCD 
sum rule method [1,2]. Clearly, an important ingredient of the QCD sum rule method is to 
determine the various vacuum condensates as precisely as possible. In the numerical pre- 
dictions of the QCD sum rule method, the largest uncertainties reside in the actual values 
for the four-quark condensates. They remained a matter of constant debate over the last 
decades [3-11]. One central point in these discussions is the question whether a four-quark 
condensate can be factorized more or less accurately into a product of two-quark conden- 
sates. The reason for this is that at present there does not exist a general approach for 
calculating this quantity. In this paper we try to present such a general approach and use 
it to analyse the factorization problem of the four-quark condensate. 

In order to make this paper self-contained, let us first recall the definition of vacuum 
condensates in the QCD sum rule formalism. In the QCD sum rule, one often postulates 
that quark propagators are modified by the long-range confinement part of QCD; but the 
modification is soft in the sense that at short distances the difference between the exact and 
perturbative propagators vanishes. To formalize this statement, one can write the "exact" 
propagator G(x) as the vacuum expectation value of the following T-product in the "exact" 
QCD vacuum |0): 

G y (x,y) = (0|T[ ft ( a ;)^(y)]|0>. (1) 
According to the spirit of Wick theorem, one formally writes the T-product as the sum 

T[qi(x)qj(y)] = qi{x)qj{y) + : q^q^y) : (2) 

S v ' 

of the "pairing" and the "normal" product. The "pairing" is just the expectation value of 
the T-product over the perturbative QCD vacuum |0) 

qi(x)qAy) = (0|r[ ft (x)«,-(y)]|0> = G^\x,y), (3) 

S v ' 

i.e., the perturbative propagator. Here we emphasize that Wick theorem is obtained in free 
quantum field theories and there is little reason to expect that it is still of validity in the 
case of exact QCD vacuum (since the vacuum |0) is highly nontrivial, it is not clear how 
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to define the normal product of field operators in this case), therefore Eq. (2) should be 
understood as a formal definition of the "normal" product. By this definition, we have the 
following two-quark vacuum condensate (0| : q(x)q(y) : |0): 

(0| : q t (x)q 3 (y) : |0> = (0|T[ft(x)^-(y)]|0> - (0|T[ft(aO(&(y)]|0) 

= (-) [Gjiiy, x) - G^y, x)} = (-) J -^e iq<y - x) [G(q 2 ) - G^( q 2 )} .. . (4) 

Here we would like to stress again that the bracketting colons : : is only a notation which 
means that we subtract the contribution of the perturbative term G pert (x) from G(x) [12,13]. 
Thus, our assumption (0| : qq : |0) ^ is equivalent to the statement G(x) ^ G pert (x). 
Similarly, we have the four-quark vacuum condensate: 

(0| : q(x)A^q(x)q(y)A^q(y) : |6) 

= (0|T [q(x)A^q(x)q(y)A^q(y)} |0> - (0|T [q(x)A^q(x)q(y)A^q(y)} |0> 

-<0| : q(x)A^ q(x)q(y) A< 2 >g(y) : |0> - (0| : q(x)A^q(x)q(y)A^q(y) : |0) (5) 

S v ' " v ' 

-(0| : q(x)A^q(x)q(y)A^q(y) : |6> - (0| : q(x)A^q(x) q(y)A^q(y) : |6>, 

V v ' S v ' 

here the stands for a matrix in Dirac and color space. 
To be more concrete, in the following let us take the case 

A (i) = A (2) 

= 7m ® Ic (where I c 

denotes the unit matrix in the color space) as an example to illustrate the general approach 
for calculating the four-quark condensate. When A* 1 ) = A( 2 ) = 7^ <g) I c , it can be shown 
that the last four terms on the right hand side of Eq. (5) vanish. For instance, for the third 
term, one can write 

(0| : q(x)A^q(x)q(y)A^q( y ) : |0> = (0| : qi(x)A\f qj(x)q m (y) A$ n q n (y) : |0> 

V V 

= ^G]Z\x,y)A^\ : q z (x)q n (y) : |0> = ^Gf^y)A^M^\ : q(x)q(y) : |0> 
= ^Tr yC [AWG^\x,y)A^](b\ : q(x)q(y) : |0>, (6) 

where the trace operation is on the color and Dirac indices. For A^ = A^ 2 ^ = 7^ <g> I c , one 
has 

7> 7C [ 7m ® I c G per % ® I c ] = ATr jC [G pert ] = 0, 

where we have assumed that we are working in the chiral limit. It can be shown that, in the 
chiral limit, the scalar self-energy function of the perturbative quark propagator vanishes to 



3 



all orders of perturbation theory [14], hence this term vanishes. For the fifth term, one can 
write 

(0\:q(x)AWq(x)q(y)AWq(y):\0) 

S v ' 

= (0| : qi(x)A^ qj (x)q m (y)^lqn(y) : |0> = -G*?*(x, x)Ag } A^(0| : q m (y)q n (y) : |0> 

V v ' 

= -^Tr 7C [G^ t (x,x)A( 1 )]Tr 7C [A( 2 )](0| : ^y) : |6>. (7) 

For A^ = A^ 2 ) = 7 M ® one has Tr 7 c[A^ 2 )] = 0, therefore this term vanishes. By similar 
procedures, one can prove that the other two terms also vanish. Therefore one has 

(0| : q(x)-f ll q(x)q(y)-f ll q(y) : |6) 

= (0\T[q(x)^q(x)q(y) lliq (y)] |0> - (0|T [q(x)^q(x)q(y)^q(y)] |0). (8) 

As can be seen from Eq. (8), in order to calculate the four-quark condensate, one must 

know how to calculate (0|T [q(x)^q(x)q(y)j^q(y)] |0) and (0|T [q(x)j fJ ,q(x)q(y)j IJ ,q(y)] |0) in 

advance. The main task now is to determine the vacuum expectation value (VEV) of the 

above T-product of four quark field operators in a consistent way. In the following we shall 

derive a general formula for these VEV using the background field method. 

Let us now study the linear response of the dressed quark propagator in the presence of 

the variable background field. The presence of the variable vector background field implies 

that the dressed quark propagator <?[V](x) (in the presence of the background vector field 

V M (x)) is evaluated with an additional term AS = / d A x q{x)i r f tl q{x)V tl {x) added to the 

usual QCD action. In the Euclidean space and the chiral limit, it can be written as 

JVqVqVA gj(a;)gj(0) exp {-S[q, q, A] + J d 4 x g(x)ry ^q(x)V '^(x)} 
yij[Vl[X) fVqVqVAe W {-S[q,q,A}+fd*xq(x)i^q(x)V„(x)} ' 1 ] 

where 

S[q, q,A]=J d 4 x ig^4 a )q + ^/;} , (10) 

and F^ v = d^A® — d p A a ^ + gf abc A^A^. We leave the gauge fixing term, the ghost field 
term and its integration measure to be understood. In this paper we adopt the following 
conventions: The Dirac matrices are Hermitian and satisfy the algebra [7^,7^] = 25^ and 
7 5 is defined as 75 = 71727374 so that Tr 7 [757^7^] = 4e^ p(T , ei 23 4 = 1- 
If one focuses only on the linear response term of C?[V](x), we have 

JVqVqVA &(x)^(0) [1 + / d 4 y q{y)i-y ^q(y)V ' M (y) + • • •] exp{-S[q, q, A]} 



JVqDqVA [1 + / d^yqiy^q^V^y) + • • •] exp{-%, q, A]} 
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1 + / ^ 2 /(0|T[g( 2/ )i 7 ^( 2/ )V / ,( 2/ )]|0) + • • • 
= (0|T[ ft (x)fc(0)]|0) + /A(6|Tfe(x)g j (0)g(y)^( 2 /)V M ( 2 /)]|0) 

-<0|T[ ft (x)$(0)]|0) J d*y{d\T[q(y)i llx q(y)V,(ym + •••■ (") 

The third term on the right hand side of the last line of Eq. (11) vanishes since 
(0\T[q(y)i^^q(y)]\0) = (since the vacuum state is Lorentz invariant, the vacuum expec- 
tation value of any local operator which carries a four- vector index vanishes identically), 

so 

Ql(x) = J A(0|Tfe(x)g j (0)?( 2 /)i 7 ^(2/)V M (y)]|0) (12) 

is the linear response term of the exact quark propagator in the presence of the vector 
background field. At this point, one remark is needed. The above background field method is 
general and one can study the linear response of the dressed quark propagator in the presence 
of background field of any other type. In all cases, up to first order of the background field, 
the expression for the dressed quark propagator contains three terms, similar to the right 
hand side of the last line of Eq. (11). In the case of the tensor, axial vector or pseudoscalar 
background field, using similar arguments as the vector case, one can show that the linear 
response term of the dressed quark propagator is given by an expression similar to (12). 
However, in the case of the scalar background field, the (0\T[q(y)q(y)}\0) does not vanish 
and the linear response term of the dressed quark propagator consists of two terms: 

= -/ d'y(mq i (x)q j (0)q(y)q(y)S(y)]\0) 

+(0\T[q i (x)q j (0)]\0) J d 4 y(0\T[q(y)q(y)S(y)]\0), (13) 

where S(y) denotes the variable scalar background field. As will be shown below, the 
contribution of the second term on the right hand side of Eq. (13) should be taken into 
account in the calculation of the four-quark condensate (0| : q(0)q(0)q(0)q(0) : |0). 

Now let us turn back to the discussion of the vector background field case. In order to 
obtain a exact expression for / d 4 y(0\T[q(x)q(0)q(y)i'j IM q(y)V^(y)}\0), we expand the dressed 
quark propagator £/ _1 [V] in powers of V as follows 

Q- l [V] = g- l [V] |v M =o +^y^ k=o V + • • • = G- 1 + V,T, + • • ., (14) 
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which leads to the following formal expansion 



Q[V] =G-GV^T fl G + 



(15) 



with 



r^(y 1 ,y 2 ;z) 



5G- 1 [V](y 1 ,y 2 ) 



(16) 



v M =o 



In coordinate space the dressed vector vertex T^(x, y; z) is given as the functional derivative 
of the inverse quark propagator £? _1 [V] with respect to the variable background field V M (x). 
Note that Eq. (15) holds for both the exact quark propagator and the perturbative one. 
Eq. (15) is a compact notation and its explicit form reads 



Q[V](yi,y 2 ) = G(y 1 ,y 2 )- J d 4 u x d 4 u 2 d 4 z G(y 1 ,ui)T^(u 1 ,u 2 ; z) ■ V^(z)G(u 2 ,y 2 ) 
= G( yi ,y2)-JdzJ—-J 



+ 



(2tt) 4 7 (2tt) z 



e -i{P+%)-yi e i{P-l)-y2 e iq-z 



*G{P+\)T,{P-q)-V»(z)G(p- q -) + 



(17) 



where P and q are the relative and total momentum of the quark-antiquark pair, respectively. 

Comparing Eq. (11) with (17), we obtain the linear response term of the exact quark 
propagator Q v (x) in the presence of variable background field 



Qlix) = J d 4 y(0|r[ ft (x)^-(0)«(y)<7M«(y)V M (y)]|0> 



(18) 



d 4 z 



d 4 q 



d 4 P 



-i(P+$)-x e iq-z 



G{P+ q -)Y,{P-q)-V,{z)G{P- q -) 



(2tt) 4 7 (2^) 4 

Setting x = and taking V M (y) = V^5(y) (V^ is an constant vector external field) in Eq. 
(18), and then multiplying by (7^)^, we have 



(6|T[g(0)7 M g(0)g(0) 7M g(0)]|0) 

r d 4 q r d 4 P 
~ J (2^)4 J (27)" 4 



Tr, 



7MG(P+|)r M (P;9)G(P-|) 



(19) 



Similarly, we have 



(0|T[g(0) 7 ^(0)g(0) 7 ^(0)]|0) 

r d 4 q r d 4 P 
= " J J^Y 1 7o^4 



Tr. 



7 ,G pert (P + |)r^ eri (P; q)G pert (P - §) 



(2^) 4 7 (2tt) 4 

where T^ er *(P; q) is the corresponding perturbative vector vertex. 



(20) 
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Putting Eqs. (19) and (20) into Eq. (8), one obtains the main result of this paper, namely 
the general formula for the four-quark condensate 



(0| : g(0)7 M g(0)g(0)7 M g(0) : |6> 




/ (2^ Tr " c r G{p + 2 )rv(F; q)G(p ~ 2 } 

\ Tr lC [ 7 ,G pert (F + |)r^ ert (P; q)G^\P - |) 




d A q f d A P 



(21) 



At this point, it is interesting to compare Eq. (21) with the general formula for the vector 
vacuum susceptibility x V [15]: 



X V oc / |^ Tr lC {^G(P)F,(P; O)G(P) - 1 ,G^\P)Y^\P ] 0)G pert (P)} . (22) 



In the QCD sum rule external field method, one often introduces various vacuum suscepti- 
bilities to characterize the properties of the QCD vacuum. It is generally believed that the 
vacuum susceptibilities and the vacuum condensates are independent quantities. However, 
comparing Eqs. (21) and (22), it is easy to see that there exists some intrinsic relations 
between the two quantities (0| : <?(0)7 M g(0)g(0)7 M g(0) : |0) and \ V ■ 

The whole procedures above are presented in the unrenormalized languages. Once the 
technique for deriving the four-quark condensate is known one can easily repeat the whole 
derivation in the renormalized languages. Essentially one only needs to modify the action 
of the theory to include the necessary counterterms and repeat the above procedure (more 
details can be found in Refs. [16,17]). 

Now let us turn back to formula (21). We note that formula (21) is formally model- 
independent, and due to presence of the subtraction term (the second term on the right 
hand side of Eq. (21)) the four-quark condensate is free of ultraviolet divergences, since 
in the large momentum region G(p) and T^k; P) coincide with G pert (p) and T^ 1 (k; P) , 
respectively. However, in order to apply this formula, one should know the values of the 
quark propagators G(p), G pert (p) and the vertex functions T^(k; P), Y p ^ rt (k; P) for the whole 
momentum region in advance. At present it is very difficult to calculate the above functions 
for the whole momentum region from first principles of QCD. Thus in practical calculation 
of the four-quark condensates one has to resort to nonperturbative QCD model. Over the 
past few years, considerable progress has been made in the framework of the rainbow-ladder 
approximation of the Dyson-Schwinger (DS) approach [16-18], which provides a successful 
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description of various nonperturbative aspects of strong interaction physics. Just as was 
pointed out in Refs. [15,19,20], the DS approach provides a general framework for calculating 
the four functions G(p), G pert (p), T^(k; P) and T p ^ rt (k\ P). For details we refer the readers 
to the above references. 

In the above we have derived the general formula (21) for the four-quark condensate 
(0| : 9(0)7^9(0)9(0)7^9(0) : |0). One can use the same procedure to derive the general 
formula for the four-quark condensate of any other type. For example, for the four-quark 
condensate (0| : 9(0)9(0)9(0)9(0) : |0), we have the following formula 



(0| : q(0)q(0)q(0)q(0) : |6> - [<0| : 9(0)9(0) : |6>] 2 
d A g r d A P 
(2^7 (2^" 4 



-i) i I I ^ Tr lC 



G(p + |)r(p ;9 )G(p-|; 



/ 70* / T ^ c [ GPeH{p + |) rpert ( p ; i) GPert ( p - f ) 



(23) 



(2-kYJ (2tt) z 

where T(P; q) and Y pert (P; 9) denote the corresponding exact and perturbative dressed scalar 
vertex, respectively. As was noted before, the linear response term Gfj{x) consists of two 
terms (see Eq. (13)), and the contribution of the second term in Eq. (13) leads to the 
appearance of the extra term — [(0| : 9(0)9(0) : |0)] 2 on the left hand side of Eq. (23). 

As was said above, one central point in the discussion on the four-quark condensate is 
the question whether a four-quark condensate can be factorized more or less accurately into 
a product of two-quark condensates. Now we have obtained a closed formula for a general 
type of four-quark condensate, we can make use of it to analyse the factorization problem. 
If one adopts the vacuum saturation assumption, one has 

(0| : qA^qqA (2) q : |6> » [Tr lC K^Tr lC K^ - Tr, c (AWA^)] [(0| : qq : |0)] 2 . (24) 

i.e., the four-quark condensate is factorized into a product of two-quark condensates. In the 
following we shall still use the four-quark condensate (0| : 9(0)7^9(0)9(0)7^9(0) : |0) as an 
illustration to analyze this problem. We note that the two-quark condensate (0| : 99 : |0) 
is related to the two-point correlation function, i.e., the quark propagator, while it can 
be seen from Eq. (21) that besides the two-point correlation function, the three-point 
correlation function, i.e., the fully dressed vertex r M , also enters into the expression for the 
four-quark condensate (0| : 9(0)7^9(0)9(0)7^9(0) : |0). Therefore, factorization of the four- 
quark condensate is essentially equivalent to neglecting the influences of the corresponding 
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dressed vertex (three-point correlation function). The simplest approximation scheme for 
the dressed vertex is the bare vertex approximation: T^(P;q) rts — ry^, T p ^ rt (P\q) — 27^. 
Interestingly, it can be shown that only under this approximation does the factorization of 
the four-quark condensate hold. The following is the proof of this assertion. 

In the case of bare vertex approximation, the first term on the right hand side of Eq. 
(21) can be written as 

d 4 q f d 4 P 



H /(2^/(2^ T ^ 
r d 4 k r dH 



(2tt)< 
d 4 k 



(27T 



I Tr ic 



[7„G(fc) 7/ xG(0] 



r drk r 
J J2^J 



dH Tr yC {-y^-i-y • kA(k 2 ) + B(k 2 )]^[-t 7 ■ IA(1 2 ) + B(l 2 )]} 



= -48 
1 

" ~3 



(2vr) 4 

d 4 k 



(27T) 



B(k 2 



[A 2 (k 2 )k 2 + B 2 {k 2 )][A 2 (l 2 )l 2 + B 2 (l 2 )} 

1 d 4 k f dH 24k ■ lA{k 2 )A{l 2 



(2ir) 4 A 2 (k 2 )k 2 + B 2 (k 2 



d A k 
W) 4 



Tr lC [G{k)\ 



(2tt) 4 J (2tt) 4 [A 2 (k 2 )k 2 + B 2 (k 2 )}[A 2 (l 2 )l 2 + B 2 (l 2 )} 

(25) 



1 f(6|:gg:|0> 12 



where we have used the fact that the integrand of the second term on the fourth line is an 
odd function of the momenta and hence this term vanishes. In addition, we have 

d 4 q f d 4 P 



(2n) 4 J (2vr 



Tr. 



^ G pert {p+ ± ) ^ GP ert (p _± ] 



0. 



(26) 



Putting Eqs. (25) and (26) into Eq. (21), it is easy to find that vacuum saturation as- 
sumption, i.e., factorization of the four-quark condensate holds in the case of bare vertex 
approximation. 

Similarly, for the case of the four-quark condensate (0| : q(0)q(0)q(0)q(0) : |0), if one 
adopts the bare vertex approximation T(P; q) m I, Y pert (P] q) ss I in Eq. (23), one can also 
show that factorization of (0| : q(0)q(0)q(0)q(0) : |0) holds. For the case of other types of 
four-quark condensate: (0| : g(0)a^g(0)g(0)a^g(0) : |0), (0| : q(0h^q(0)q(0)l^q(0) : |0) 
and (0| : ^(0)755(0)5(0)75^(0) : |0), one can also prove that factorization holds in the case 
of bare vertex approximation. 

From the above results, we reach a general conclusion: by means of the background field 
method, one can derive general formula for all five types of four-quark condensates, which 
contains the corresponding type of fully dressed vertex, and factorization of the four-quark 
condensate holds only when the dressed vertex is taken to be the bare one. 
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To summarize, by differentiating the dressed quark propagator with respect to a variable 
background field, the linear response of the dressed quark propagator in the presence of the 
background field can be obtained. From this general method, using the vector background 
field as an illustration, we derive a general formula (21) for the four-quark condensate (0| : 
g(0)7 M g(0)g(0)7 M g(0) : |0), which contains the corresponding dressed vertex T^q; P) and 
T^ rt (q; P). In the case of bare vertex approximation, it is shown that factorization of the 
four-quark condensate holds exactly. Finally, we want to stress that the background field 
approach adopted in this paper is general enough and can be applied to the calculation of 
the vacuum expectation value of the T-product of any four quark field operators. Especially, 
it can be used to study the two-point correlation functions of the mesonic current, which 
play an important roles in nonperturbative QCD and hadronic physics [21]. 
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